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The time-dependent CP asymmetry in exclusive B'^{t) K'^'j decays has been proposed as 
a probe of new physics in B decays. Recently, this method was extended to radiative decays into 
multibody hadronic final states such as B^{t) Kst^^^ and B^{t) ty^tv'j. The CP asymmetry in 
these decays vanishes to the extent that the photon is completely polarized. In the Standard Model, 
the photon emitted in 6 — > 57 has high left-handed polarization, but right-handed contamination 
enters already at leading order in A/m;, even for vanishing light quark masses. We compute here 
the magnitude of this effect and the time dependent CP asymmetry parameter Sj^^^Oj. We find 
that the Standard Model can easily accomodate values of S as large as 10%, but a precise value 
cannot be obtained at present because of strong interactions uncertainties. 

PACS numbers: 12.39.Fe, 13.60.-r 



I. INTRODUCTION 

The standard model (SM) predicts that photons are 
predominantly left-handed in b —^ (q — s,d) decay 
(and right-handed in 6 — *■ qj). In the presence of new 
physics this prediction can be changed, and a significant 
right-handed photon amplitude can appear in 6 — > 37 
decays. Several methods have been suggested for testing 
this prediction in radiative B decays [ij, I3] • 

One of these methods makes use of time-dependent 
CP violation in B'^ with / a CP eigenstate [ij. 

Since 7^ and jn cannot interfere, the time dependent 
CP asymmetry 



r[^°(0^/7]-r[i3°(t)^/7] 
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— Sfj sin( Am t) ~ C cos(Ato t) , 
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is sensitive to the ratio of the right- /left-handed photon 
amplitudes. These can be parameterized as 



rfe^'f =Vcp{f) 



(2) 



where ijcpif) is the CP eigenvalue of the state /. This 
method has been extended in Ref. 0| also to decays into 
multibody final states, such as for example B — > Kstt'^j. 

Summing over the unobserved photon polarization, the 
CP violating parameter in Eq. Q is given by 



2rf 



1 



cos Sf sin 2/3 . 



(3) 



In the SM, it is usually assumed (incorrectly, see 0) that 
fq ~ rnq/nii,, which leads to a small CP asymmetry S ~ 
2% in the b ^ transition. We used in this estimate 
sin 2/3 = 0.685 ±0.032 0. 

Measurements of the CP asymmetries in S — > K*j 
were reported by BABAR j3] and BELLE 0, see Ta- 
ble I. In addition, these two collaborations measured the 



CP asymmetry in the nonresonant mode Sj^sTrO-y in two 
different ways: a) BABAR excludes the K* resonance 
by integrating over the Kstt'^ invariant mass range 1.1 
GeV < Mk^^o < 1.8 GeV; b) BELLE includes both 
resonant and nonresonant modes by integrating over the 
range niK + my,- < Mk^^o < 1.8 GeV. The error in these 
determinations is still too large to allow a meaningful 
comparison with the SM prediction. With a view to im- 
proving the statistics of such measurements, we would 
like to assess the feasibility of combining the resonant and 
nonresonant measurements. Also, it is clear that search- 
ing for new physics with such measurements requires a 
reliable estimate of the standard model background. 

At leading order in l/mt, the photon emitted in _B — > 
fC*7 is always left-handed polarized, to all orders in as 
P|. However, in multibody decays such as B ^ Kstt'^j a 
right-handed component appears, already at leading or- 
der in l/r7ifc, in the kinematical region with an energetic 
kaon and a soft pion. This is mediated by a B* pole 
diagram, with the B* — > K^j^ amplitude calculable at 
leading order in A/rUb using factorization in QCD as fol- 
lowing from Soft CoUinear Effective Theory (SCET) ^. 
The presence of two hadrons in the final state evades the 
helicity argument which forbids a right-handed photon in 
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-0.21 ± 0.40 ± 0.05 +0.01 ± 0.52 ± 0.11 
0.9 ±1.0 ±0.2 +0.20 ±0.66 
-0.06 ±0.37 +0.08 + 0.41 + 0.10 



TABLE I: Experimental results for the CP asymmetry pa- 
rameter Sf for / = K*^'y and / = Ks'J^'^'y from BABAR and 
BELLE. The BABAR nonresonant region includes all states 
with 1.1 GeV < rnKs-^a < 1.8 GeV, while BELLE uses the 
range 1.0 GeV < mj^^^o < 1.8 GeV. The errors shown are 
statistical and systematic, respectively. 
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B K*j. Note that a right-handedphoton is present in 
inclusive B -—^ Xgj at leading order [9j. The right-handed 
photon couples to the charm quark loop induced by the 
4-quark operator O2 = (sc)(c6), which gives equal rates 
for b sgjL and b — > sg^R. In addition to this leading 
order effect, a significant right-handed photon amplitude 
in B ^ K*"f can appear at subleading order in A/ mi, 
from graphs with photon emission from the charm quark 
loop. 

The purpose of this paper is to study in more detail the 
magnitude of the leading order effects described above on 
the time-dependent CP asymmetry in B — > KsTr'^-f- In 
Sec. II we introduce the effective theory formalism used 
in our computation. This is a combination of the soft- 
collinear effective theory (SCET) with the chiral pertur- 
bation theory recently proposed in Ref. [s^. In Sec. Ill 
the helicity amplitudes are written down, and used to 
compute decay distributions for B Kstt^j with a 
right-handed photon. Sec. IV gives the results for the 
time-dependent CP violation SKs-n-o-y in the kincmatical 
region with an energetic kaon and a soft pion. This has 
a significant overlap with the region used in the BELLE 
and BABAR measurements of the time-dependent CP 
asymmetry into a nonresonant Kstt'^ state. Sec. V 
summarizes our results. Readers interested in the phe- 
nomenology of the decay can skip the formalism and pro- 
ceed directly to Sec. III. 



II. EFFECTIVE THEORY FORMALISM 

The exclusive radiative decays B AT* 7 can be de- 
scribed in the large recoil region in factorization. At 
leading order in A/irib with A ^ 500 MeV, the existence 
of such factorization relations has been demonstrated in 
0,0,0,0,01 at lowest order in as, and proven to 
all orders in abusing the soft-coUinear effective theory 

The b transitions with an energetic s quark are 

mediated in SCETj by the effective Lagrangian 



mbc{uj) Sn,u-^'^PL bv 



(4) 



with Nq 



— — ij^,. The relevant modes are soft 
quarks and gluons with momenta kg A and coUinear 
quarks and gluons along n j29l |. We use everywhere the 
SCET notations in Ref. |2^ . We choose the photon mo- 
mentum to move along the —63 direction = E, 



such that the hadronic system has a large momentum 
along the opposite direction n^. The hard scale in this 
problem is Q = n • px^ ^ "m-b, and the expansion pa- 
rameter in SCET is ~ A/Q. We denoted n^,n^ unit 



transverse photon field is Aj^, and its polarization vec- 
tors are e+ = {^1 (right-handed photon) and 

e_ = 71 : 7f ;0) (left-handed photon). 

We neglect here and in the following s quark mass ef- 
fects, which can be included straightforwardly |23|. The 
photon coupling to the spectator quark in the B in- 
troduces new factorizable operators containing collinear 
modes along the photon momentum |29|, |30|, |31| . These 
spectator effects do not contribute to the right-handed 
photon amplitude, and we return to them below (see the 
discussion around Eq. H14|l 'l. 

The first operator in Eq. Q scales like 0(A°) and cou- 
ples only to left-handed photons. The 0(A) operators 
QiiL,R) couple to left- and right-handed photons, respec- 
tively, and are defined as 



0(i«)(c^i,c^2) = 
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n ■ V 



■V 

i0n 



PRb-a, 



4^PRb,. (5) 



At lowest order in matching, the Wilson coefficients can 
be extracted from the computations of |l3.ll^ll^ll5lll9| 
and are given by [29| 
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&lfl(t^l, W2) 
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(6) 
(7) 



4 ^^r' ) +g[g3-6,8,«.(mfc)] 
V mi / 



4 4 11 
C7 - -C3 - -C4 + -C5 + -Cq (8) 



in the operator basis of Ref. |13| ■ Beyond tree level the 
Wilson coefficients c, 6il and 61 b receive hard corrections 
~ as{mi,) from charm loops [iSj and from matching onto 
SCETj. The 0(a.^(mh)) matching corrections are known 
only for c{uj) [21|, so for consistency we do not include 
them in any of the coefficients. The function k{x) appears 
in the 3-point function with a charm loop, and is given 
by J^] 




x < 4 

(9) 

f ) x>A 



In the Wilson coefficients of the 0(A) operators 6il and 
biR we neglect small contributions from the penguin op- 
erators Os-e and the gluon dipole operator Og. 

After matching onto SCETn [i^ , the b sj effective 
Lagrangian Q contains both factorizable and nonfactor- 
izable operators 



light-cone vectors satisfying n 



0,n • n = 2. The 



(Of + 0*^^')A'' + 



(10) 
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where the elhpses stand for higher dimension operators. 
The details of this matching are given in Refs. |25ll2g . l27l |. 
and we give here only the points essential in the following. 

The nonfactorizable operators couple only to the left- 
handed photon field, e^^'OJJ^ = 0, while the factorizable 
operators couple to both left- and right-handed photons. 
Working at tree level in matching SCETj SCETn, the 
factorizable operators are 



1 

2uj 



2uo 



dzda;dfc_|_6i/j(z) J|| (a;, z, fc+) 

,CJ2 



(11) 



where we used a momentum space notation for the soft 
nonlocal operators 



dX 



r(An/2)y„(A,0)6i(0). (12) 



The functions biL{z) and biji{z) appearing here are re- 
lated to the Wilson coefficients in Eq. as bi{z) — 
— z)uj,zLj). The momentum labels of the coUinear 
fields are parameterized as loi — xu},u}2 = —(1 — x)lo. 
We denoted here with Jj^ y jet functions defined as in 
Ref. [s^. They have perturbative expansions in ai(/ic) 
with /x^ ~ QK. At leading order they are given by 



J\\{x, z, k+) = J±{x, z, fc+): 



Nr. 



xkj^ 



■5{x~ zlli) 



The 0{al{^c)) corrections to the jet fmictions have been 
recently obtained in Ref. [isf . 

Another class of factorizable operators not present in 
Eq. arise from the photon coupling to the spectator 
quarks (2^ 

(Photon coupling to the final state 
quarks leads to power suppressed operators H^].) After 
matching onto SCETn, they are given at leading order 
in as{mh) by 



01' - 1^' E ".M) 



(14) 



q—u,d.s 



X y dk^Jsp{k^)eq{qk_jf^yiiPLbv){sn,uji-^PLq7i,uj2) 

with bspU^^) = VubV:,{Ci + C2/N,)5,u ~ V^t6K*(C4 

C3/iVc) + 0(a.(Q)) and 



to left-handed photons. For consistency with the other 
factorizable operators included, we work to 0(a°((5)), 
but keep terms of 0{as{fJ.c)) in the factorized ampli- 
tude. The dominant term ~ C1.2 contributes only to 
5° K-TT+"f, but not to 50 KsTr°-/. In Eqs. © we 
neglected the contributions from Os-e.s to the Wilson 
coefficients c(aj), biL,iB.{uJi), so for consistency we neglect 
such terms also in Eq. fTljl . 

The SCET formalism introduced above has been 
applied to prove factorization relations for exclusive 
semileptonic B M£i/ and radiative B M7, M£~^£~ 
decays into one energetic light hadron, with M a ligh t 
pseudoscalar or vector meson [j I2I Ei. EtI lil liil Isoj . 
In all these cases, the transition matrix element is written 
as the sum of a soft (nonfactorizable) and hard-scattering 
(factorizable) terms, as follows. 

The matrix elements of the nonfactorizable operators 
are parameterized in terms of soft form factors. In our 
calculation we require only the matrix element 



{K*ip, r,)\sJ*^PLb,\B{v)) = {e* 



(16) 



where we use the SCET/ notation for the operators ob- 
tained from them by matching onto SCET//. 

The matrix elements of the factorizable operators O^^'^* 
given in Eq. (|lll) are given by convolutions of soft and 
coUinear matrix elements with the Wilson coefficients. 
The matrix elements are parameterized in terms of light- 
cone wave functions of the B and K^*^ mesons |2^. In 
the B K* transition, only the left-handed photon am- 
plitude is nonvanishing 9, 29], and is given by 

^fact(^ ^ K*^^^ = {K*jL\n.s\B) (17) 
= NomtmB / dz6iL(z)Cff '(z) 

JQ 



where the factorizable function (f]^ (z) is defined as 



ruB 



dxdk^J±(x, z, A;_|_)(/):^(fc+)(/);^. (x) 



(18) 



In Ref. [38j it was proposed to extend the application of 
this formalism also to multibody B decays to final states 
containing one energetic meson and one soft hadron. We 
summarize briefly the main points of this approach, be- 
fore proceeding with the details of the computation. 

The matrix element of the nonfactorizable operator is 
parameterized in a manner similar to Eq. Hlt)|) in terms 
of a new soft function 

H'1\B ^ K^7rjL)^Nom,cin-pK^)Cf'''' {Ek, E^) (19) 



Jsp {k— ) 



k^ 
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with L = 1ok[(- 
to 0{as{nc)) m 



2E^k 



An 



1 



y)) (15) 



ie) I /i^] a jet function known 
The spectator operator couples only 



This nonfactorizable amplitude couples only to left- 
handed photons, just as in the case of the B K*^ 
transition. Furthermore, the soft function Q^^"^ is re- 
lated by a symmetry relation to a similar function ap- 
pearing in multibody semileptonic decay B — > TTTrtv [38| , 
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analogous to the appearance of a common nonfactoriz- 
able amplitude Cf^^ in both rare and semileptonic form 
factors llSIIllii. 

The matrix elements of the factorizable operators in 
Eq. are also given by convolutions of hard, jet and 
soft factors, as in the case of the B K* transition dis- 
cussed above. At leading order in A/Q, soft and coUinear 
modes decouple in the SCET Lagrangian [i^, which is 
the statement of soft-collinear factorization. This fact 
has two important implications. First, the soft pion does 
not couple to the coUinear meson in the final state M . 
Second, the matrix elements of factorizable operators cor- 
responding to the transition B — > M„7r factor as 



(i^„7r|Ofact|S) = (X„|Oc|0) X {^\Os\B) 



(20) 



and are calculable in terms of the kaon light-cone wave 
functions and a new soft matrix element of the Os oper- 
ator in the B —f tt transition. 

The soft operator Os required here appears in the biji 
factorizable operator. We define it as 



dX 

£ 

4tt 



+ g(A^)y„(A^,0)^+^Pfl6,(0)(21 



Its i? — > TT matrix element is parameterized in terms of 
one soft function S{k+^t^ defined as 

{AP^)\Os(k+)\B{v)) = -2(e+ ■p^)S{k+,t\0 (22) 

with t — niBV — Ptt and ^ = n -Ptt/ {n ■ v). The support of 
this function is —n ■ Ptt < k^ < oo. This is the analog for 
B physics of the generalized parton distributions (GPD), 
commonly encountered in nucleon physics. 

The complete factorization relation for the right- 
handed photon amplitude can be now written down as 

H+{B ^ Kttjb) - 2iVo/A'(e+ • P.) / dzdxbiR{z) 



(23) 



X / dfc+ J|| (x, z, A;+)5(fc+, i X)4>k{x) 



The predictive power of such relations depends on the 
existence of reliable information about the soft function 
S. Eventually the function S should be extracted using 
B decays data, or constrained by lattice QCD compu- 
tations. In the soft pion region, the soft function S is 
fixed by chiral symmetry in terms of one of the B meson 
light cone wave functions [s^- We will use in this pa- 
per the result for S predicted at leading order in chiral 
perturbation theory. 

The predictions of chiral symmetry for the couplings 
of Goldstone bosons are most conveniently derived us- 
ing chiral perturbation theory. The applicability of this 
approach is limited to problems describing only soft 
hadrons. The extension to heavy hadrons is possible, pro- 
vided that the large scale rub is eliminated by going over 
to HQET. The corresponding chiral effective theory is 



the heavy hadron chiral perturbation theory (HHChPT) , 
and its degrees of freedom are heavy meson s pin doublets 
H = {B,B*) and the Goldstone bosons [Ulilil. 

The effective Lagrangian that describes the strong in- 
teractions of the Goldstone bosons with the ground state 
heavy mesons is 

mills 



c 



(24) 



+ izgTr [H^Hb^^js] [^^5^^ - ^d'^^^] + • • • 

where the ellipsis denote light quark mass terms, 
0(l/mb) operators associated with the breaking of heavy 
quark spin symmetry, and terms of higher order in the 
derivative expansion. The pseudo- Goldstone bosons ap- 



pear in the Lagrangian through ^ 
where 



n 



TT 

K- 



V6 



K+ \ 



e) 



(25) 



75 V 



with the pion decay constant /jr — 135 MeV. These fields 
transform as 



(26) 



under chiral SU{3)l x SU{3)ii transformations. The su- 
perfield Ha contains the pseudoscalar and vector heavy 
meson fields Ba and B*^ with velocity label 



K^7^ - Bal5] ■ 



(27) 



The flavor index runs over a = 1,2,3 corresponding to 
Ba = {B-,B",Bs). Under chiral ^[/(S)^ x SU{3)r, the 
superfield Ha transforms as 



Ha Hb Uf 



ha ' 



(28) 



The numerical value of the coupling g = 0.5±0.1 is taken 
to cover a range compatible with its determination from 
D* -> Dtt decays g = 0.59 ± 0.08 "io"! and lattice QCD 
.g = 0.48 ±0.03 ±0.11 0> 5 = 0-42 ± 0.04 ± 0.08 42]. 

We consider next the matrix element {■k\Os\B) in chi- 
ral perturbation theory. This requires the chiral repre- 
sentation of the nonlocal soft operators Os- Consider 
light-cone nonlocal heavy-light bilincars of the form 



oiAk, 



dx_ 
47r 



r (a:- )r„ {x- , Q)PR.,LTh, (0) . 



(29) 



Under the chiral group they transform as (3^,1^?) and 
(1l,3/j), respectively. For each case, there is a unique 
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operator in the effective theory with the correct transfor- 
mation properties jSSj 

02(fc+) = -Tt [aL{k+)Pi,TH,ill (30) 

Ol{k+) = ^Tr [aR{k+)PLTHb^^,a] (31) 

The common matrix dL(fc+) — aji{k-^) — a{kj^) is given 
by 

a{k+) = /i3V^[#?(fc+) + (32) 

where 0^(fc+) are the usual hght-cone B meson wave 
functions, defined by 

/ ^e-*'=+-(0|g,(z_)r„(z_,0)6i(0)|B(t;))= (33) 

-^fsms <^—^[^n ■ V(j)f{k+) + ijln ■ v(j)^{k+)]"fr-, j 

The operators in Eqs. l|5n|) . with can be used to 
compute the matrix elements of Ol.b on states with a B 
meson and any number of pseudo Goldstone bosons. In 
particular, they give the following prediction for the soft 
function S{k^,P,Q defined in Eq. ()22|l at leading order 
in the chiral expansion (38| 

^(fc+,t^c) = ^— ^<^?(fc+) (34) 

We will use this result in the numerical evaluations of 
this paper. No information can be obtained using chiral 
symmetry about the S function for —n • < fc+ < 0. 
For soft pions this contribution is likely to be small, so it 
will be neglected in the following. 

Finally, we comment briefly on previous applications 
[3^ of HHChPT to B decays into multibody final states 
containing soft Goldstone bosons. These applications in- 
volve the generalization to B* decays of the factorization 
formula for nonleptonic B decays, with the B* appear- 
ing in intermediate states of pole diagrams. The usual 
HHChPT [3^ [3^ [33 methods are apphed to compute 
the pion coupling in both pole diagrams and in contact 
diagrams. 

There are several issues with such a simplified ap- 
proach: i) the application of chiral perturbation theory 
to the nonfactorizable operators Onf contributing to the 
_B — !■ TT transition with an energetic pion is problematic. 
Since these operators couple to both soft and coUinear 
modes, loop corrections to their matrix elements do not 
have a well-behaved power counting. In our approach 
these contributions are simply parameterized by new soft 
functions ^j^^'^ , which are related by symmetry relations 
to similar matrix elements appearing in other processes, 
ii) the pion contact terms, such as that in Fig. 2b, can 
be computed only for the factorizable operators, (but not 
for the entire weak vertex), and are given by factorization 
relations. 
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FIG. 1: The phase space of the B° — > Kstt^J decay in 
variables {Mk^, Et^). The vertical line denotes the maxi- 
mum Kn invariant mass used in the BABAR and BELLE 
measurements. The 3 regions shown correspond to (I) soft 
pion Et^ ~ A; the shaded region E-^ < 0.5 GeV (implying 
Ek > 2.18 GeV) shows the region of applicability of ChPT; 
(II) coUinear pion and kaon Ett ~ Q,Ek > 1 GeV; (III) soft 
kaon Ek <1 GeV (implying > 1.7 GeV). 

In the next section we will combine the pieces of the 
factorizable amplitudes, add in the nonfactorizable am- 
plitude and write down the complete result for the multi- 
body B —^ Kn^ amplitudes. 



III. HELICITY AMPLITUDES AND DECAY 
RATES 

We will use the formalism described in Sec. II to com- 
pute the amplitude for the decays B Kstt'^'J and 
B — > K~TT~^j in the kinematical region with one ener- 
getic (coUinear) kaon and one soft pion. To establish the 
region of validity of our computation, we show in Fig. 1 
the phase space for this decay, in variables {Mktt, E,^), 
with M|.^ = {pK + P-nf- 

We distinguish three distinct regions for the pion and 
kaon energies in i? — > Ktt^ decay (see Fig. 1): 

I) ~ A, Ek ^ Q 

II) E^^Q,Ek-Q 

III) E^^Q.Ek^K 

These three regions are treated differently in the 
SCET, and the heavy quark mass scaling of the decay 
amplitudes is correspondingly different in each of them, 
as follows. 

The region (I) contains a soft pion and an energetic 
kaon. Part of this region, but not all, can be treated using 
the SCET -t- ChPT combination considered in this paper. 
We subdivide it into the soft pion region with iJ^r < 0.5 
GeV, where chiral perturbation theory is valid, and the 
intermediate pion region 0.5GeV < i?7r < 1-0 — 1.5 GeV 



(which we wih call, for lack of a better name, the hard- 
soft pion region). 

The region (II) includes coUinear pion and kaons. In 
general this configuration can have the kaon and pion 
momenta moving in different directions forming a large 
angle in the B rest frame Ot^k ^ 0{1). In our case, the 
experimental constraint Mktt < 1-8 GeV forces the angle 
to be small O^k ~ 0{A/Q) (valid for ~ AQ), such 

that the tt,K constituent partons can be described by 
coUinear fields with a common n. 

The region (III) contains an energetic pion and a soft 
kaon Ek ^ A, ~ Q. This region is not described 
by the leading order SCET operators in Sec. II, and the 
corresponding amplitudes are suppressed by at least one 
power of A/Q relative to those in regions (I) and (II). 

We start by writing down the helicity amplitudes for 
the B K*j decay at leading order in A/m^, by com- 
bining the partial results in Sec. 2. At this order, only 
the left-handed photon amplitude is nonvanishing 



(35) 
(36) 



: (c(ms )Cf 



BK* 



dzb,Liz)Cff {z) 



NombmBCfgfiO) . 



We defined here the effective tensor form factor g'^{0), 
which absorbs the contributions of the operators other 
than O7. Similar factorization relations are expected to 
hold also for the B Ktt^ transition in region (II), 
with the K* light-cone wave function replaced by a two- 
body Kt: light-cone wave function. We do not pursue 
this further here, but note only that the vanishing of 
the right-handed photon amplitude at LO observed in 
Eq. should hold also for B i^„7r„ in the region 
(II). Technically, this follows from the vanishing of the 
B — + matrix element of the soft operator in Eq. 
multiplying biR. 

The amplitudes for the multibody transition B — > 
Kt:^ in region (I) are given at leading order in A/mh 
and in the chiral expansion by the graphs in Fig. 2 

(I) : H+{B ^ Kttjr) ^ No^mlSRip^) (37) 

BK / 



H^{B ^ Ktt^l) = Non-pKcimBK^''''iEK,E^) 

(38) 

The nonfactorizable operators O^j contribute only to 
the left-handed photon amplitude, and the corresponding 
matrix element is parameterized by Cj^^^. On the other 
hand, the factorizable operators Ol^^^^ contribute only to 
the right-handed photon amplitude. (We consider only 
B'^ KsTT^^ decays, for which the spectator factorizable 



B B* ^ K B B 



^ 7t 

(a) 



n 7t 

(b) (c) 



FIG. 2: Diagrams showing leading order contributions to the 
decay B —> Kn^ with one coUinear kaon and a soft pion. The 
fiUed circle in (a) and (b) represents a factorizable operator 
Ofact, while the filled square in (c) represents a nonfactoriz- 
able operator Onf. 



operators Ospoc do not contribute in the approximation 
used here of neglecting Oa-e-) 

Note the appearance of a nonvanishing right-handed 
photon amplitude at leading order in the region (I). This 
amplitude is factorizable and can be computed as ex- 
plained in Sec. II. The HHChPT diagrams required for 
its computation are shown in Fig. 2a. The factorizable 
function C^^(z) appearing here is defined in analogy 
with the function in Eq. H18|) . with the replacements 
fk' ~* Ik ,(t>K*{x) 4>k{x),Ji_ — > Jj|. The depen- 
dence on the pion momentum is contained in the soft 
functions Sr{p-„) given by 



SR{V-n) 



9 g+ ■ Ptt 
UE^ + A 



(39) 



with A = niB- - mB = 50 MeV. 

Finally, in the kinematical region (III) with one soft 
kaon and an energetic pion, the effective Lagrangian 
Eq. Q does not apply. The leading SCETj operator 
mediating such a transition contains the sTby soft cur- 
rent, with at least two insertions of the soft-coUinear La- 
grangian, acting on the spectator quark 



T{(.-r6„),z£«,z£«} 



(40) 



This is suppressed by at least ~ A/Q relative to the 
operators in Eq. which implies that the decay ampli- 
tudes in this region must be power suppressed relative to 
those in regions (I) and (II). 

We summarize the different contributions enumerated 
above by showing them in graphical form in Fig. 2. We 
emphasize our different treatment of the factorizable and 
nonfactorizable operators: the matrix elements of the 
factorizable operators are computed in chiral perturba- 
tion theory, and include the B* pole and contact terms 
(Fig. 2(a), (b)). The matrix element of the nonfactoriz- 
able operator (Fig. 2(c)) is parameterized in terms of a 
new soft function (^f^^"^ . This new soft function appears 
only in the left-handed photon amplitude, and is related 
by symmetry relations to a similar soft function which 
can be determined in principle from B 7r„7r£P [s^ . 

The only region where the right-handed photon am- 
plitude B^ — > KsTT^^ contributes at leading order in 
A/Q is the region (I) with a soft pion. In region (II) 
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the right-handed amplitude is suppressed by A/Q rela- 
tive to the left-handed amplitude, and in region (III) both 
amplitudes are suppressed by at least A/Q. The contri- 
bution of the region (I) with a hard-soft pion 0.5 GeV 
< Ej^ < 1.0— 1.5 GeV will be estimated by assuming the 
validity of HHChPT in the entire region (I). We proceed 
to compute the right-handed photon effect in region (I) 
on the decay rates and time-dependent CP asymmetry. 

The left-handed photon amplitude H^{B Ktt^l) in 
Eq. H38(l (region (I)) docs not include the K* pole contri- 
bution, although this likely dominates numerically in the 
resonant region Mk-r ~ Mk- ■ This contribution is para- 
metrically suppressed by A/m;, in the soft pion kincmat- 
ical region. The reason for this is that by soft-coUinear 
factorization at leading order in A/to^, the coupling of 
a soft pion to two coUinear hadrons, K*Knns, must be 
power suppressed. On the other hand, the K* pole con- 
tribution is numerically enhanced by the K* propagator, 
so we will include it in our computation, despite being 
formally of higher order in A/uib relative to the latter. 

In the absence of data on C±^^ j we will model it by a 
K* pole contribution. We introduce the following model 
for the —> Kttj decay amplitudes in the kinematical 
region (I) with one soft pion and an energetic kaon (this is 
similar to the model used in Ref. |4^). In the left-handed 
photon amplitude we neglect the nonresonant contribu- 
tion and keep only the K* pole term (Fig. 2c) 



H 



model 



(41) 



with H^^' = A{B K*jl) the 1-body hehcity am- 
plitude given in Eq. (ESJ, and BWif.(MK^) = [Af|.^ - 
M]^, +iMK'TK']~^ the Breit-Wigner function for a K* 
resonance. The amplitude for i?° Kstt'^Jl is given 
by Eq. (jUJ) multiplied by 1/2. The K*"K+7r- cou- 
pling with a charged pion can be extracted from the total 
K* Kit decay width, T — 5|-./f^P^/(167rm|..), which 
gives gK-KiT = 9.1. 

The right-handed photon amplitude is given by the 
sum of the K* and B* resonant terms 



H 



model / dO 



{B" ^ K~7:+^r) 



(42) 



2/^ E^ + /^ 



+ gK'Me+-p^)BWK'{MK^) 

The B'^ KsTT^jR, amplitude has an additional factor 
of 1/2. We included here also a nonvanishing resonant 
B — > K*jf( right-handed photon amplitude, which is in- 
troduced by a nonvanishing strange quark mass, and by 
power suppressed contributions neglected in Eq. H35() . We 
will parameterize it as 



H 



BK* 



nib 



(43) 



inS 1 dTji 




. 8 



1.2 



1 . 6 



(GeV) 



FIG. 3: The decay rate (xlO^) with a right-handed photon 
normalized to the B — » K*'yL rate, with a 
cut on the pion energy Et, < 0.5 GeV, and its components, 
computed as described in the text. The black curve gives 
an upper bound on the total decay rate. The dark grey curve 
shows the B* — K* interference term, and the light grey curve 
shows the contribution of the B* pole graph (magnified by a 
factor of 10"^ relative to the other curves). 



The potentially leading mechanism contributing to hs 
has been identified in Ref. 0, and it arises from charm 
loops coupling to the B and K* through soft gluons. Al- 
though no first principles calculation of this parameter 
is yet available, it can be estimated from a simple power 
counting argument as 



6 Cj nib 



(44) 



In our numerical evaluation we will use = 5%, keeping 
in mind that this estimate is on the lower side of the 
dimensional estimate. 

In the remainder of the paper we will use the model 
described by Eqs. H41l) - (|43|l to compute distributions and 
decay rates for the B Kstt'^j decay. We start by 
computing the right-handed photon rate; although not 
directly observable, this quantity will illustrate the rela- 
tive importance of the different mechanisms contributing 
to the wrong-helicity photon amplitude. 

The B Ktt^ decay rate is given by 



d^T{B Kttj) _ 1 
dE^dMl^ ~ 2(47r)3 



\H\ 



(45) 



where H± are given in Eqs. 1)41(1 and ((42|l . respectively. 
In the limit of a very narrow K* the integrations over 
(Mktt, Et^) can be performed exactly, and the well-known 
result for the B K*j rate is recovered 



dE^dMl, 



dE^dMl^ 



(46) 



E. 



(0) 



STrrrip 



\H_{B^K*-fLt 
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TABLE II: Input parameters used in the numerical compu- 
tation, and results for the effective Wilson coefficients and 
factorizable matrix elements. The values of the effective Wil- 
son coefficients are quoted at the scale fi = 4.8 GeV. The 
strange quark mass is taken from l45|. 



pole 


4.8 GeV 


C2 


1.107 




1.4 GeV 


Ct 


-0.343 


m,(2 GeV) 


78 ± 10 MeV 


fB 


200 MeV 


gfio) 


0.3 


Ik 


170 MeV 




350 MeV 


QK'Ktt 


9.1 


9 


0.5 ±0.1 


id} 

Ke ^ 


-0.013 - 0.045j 



Here iS^"-' = (m^ — m^. )/(2mB) denotes the photon 
energy corresponding to the 2-body kinematics. 

It is convenient to express the B K^tt^^h decay 
rate by normalizing it to the B K*^^ decay rate (a 
factor of 1/4 has to be added for B^ Kstt^^r) 



1 



2(47r)2£;. 



(0) 



(47) 



gne 



icf} 



4°) VAE. 



K 71 



Mi. 



with Mj^^ — — 2mBEy. The two terms give the con- 
tributions of the B* resonant pole, and that of the K* 
resonant right-handed amplitude. The hadronic dynam- 
ics in the B* resonant contribution enters through the 
RG invariant ratio 



Ke 



rriB 



mtCfgfiO) Jo 
-0.013- 0.045j 



dz6ifl,(z)Cf^"(z) 



(48) 



IV. TIME-DEPENDENT CP ASYMMETRY 



We compute in this section the mixing-induced CP vi- 
olating parameter Sxs-K^-y in the Standard Model. We 
start by defining the time-independent amplitudes 



Ar 
Al 
Ar 



H{B° ^ Kstt^jr) 
iJ(B" ^ /^stt^l) 
if(B" ^ Kstt'^jr) 



(49) 
(50) 
(51) 
(52) 



Since the 6 — > s transition is CP conserving, there are 
relations among these amplitudes, such that only two of 
them are independent. We choose the independent am- 
plitudes to he Al, Ar, and obtain the remaining two am- 
plitudes from them by CP transformations. Charge con- 
jugation exchanges particles and antiparticles, and parity 
takes 7i <-> ^r and changes the directions of momenta. 
We apply a rotation by 180° around the x axis, which 
restores the momenta to their original directions. The 
effect of the rotation is to multiply the amplitudes with 
ITiTTi (—)''' = +1 (with JJ'^ the spin-parity of the parti- 
cles), and exchange e+ ^ £_. This gives 



Al^'-±^Ar. 



Ar 



At 



(53) 



The time-dependent differential rate has a form similar 
to Eq. gSJ (with i = L,R) 



1 



2(47r) 



3 ,-(|A,P + |i.n 



dE^dMl^ 

^ ""{l + QcosAmi-S'iSinAmt} (54) 



X -e 
2 



In the numerical evaluation of this parameter we used 
the tree level result for the jet functions Eq. (|13|l and the 
lowest order matching result for biR from Eq. The 
remaining required parameters are listed in Table II. 

The result Eq. <|T7|l can be used to compute the en- 
ergy spectrum and integrated rate with a right-handed 
photon, with an upper cut-off on the pion energy. The in- 
terference of the two terms depends sensitively on the un- 
known strong phase <ps ■ For this reason, we will give only 
an upper bound on this rate, obtained by assuming that 
the two terms have the same strong phase and interfere 
constructively. The resulting photon energy spectrum 
and its components are shown in Fig. 3, for the central 
values of the parameters. For completeness, we quote 
also the fraction of events which survive a pion energy 
cut ncnt{E™^'^). This can be computed from the K* pole 
contribution to H^{B Kttj) and is: ncut(0.5 GeV) = 
11.4% , ncut(l GeV) = 50.5% , ncut(l-5 GeV) = 88.4%. 




0.8 1 1.2 1.4 1.6 1.8 



Mk. (GeV) 

FIG. 4: The time-dependent CP asymmetry parameter 
^KsTrO-y (in percent) as a function of Mk-r ■ The three lines cor- 
respond to (from bottom to top): hsCostj)^ = —0.05,0,0.05. 
We used here i;^'"" = 0.5 GeV. 
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with 



Im [e-'^'PAjA*) 
|v4.P + |A,|2 



(55) 
(56) 



From this expression, results for the CP violating coef- 
ficients integrated over parts of the phase space can be 
straightforwardly obtained. 

The BELLE and BABAR Collaborations measured the 
S and C parameters integrated over all E-„ and a range of 
Mktt- We compute the SM values of these parameters, by 
integrating the time-dependent distribution Eq. 154(1 with 
appropriate cuts. Applicability of the chiral perturbation 
theory computation of the Af( amplitude requires that 
we restrict the pion energy in the B rest frame by E^^ < 
with E'^^'''' = 500 MeV. 

Integrating over E^^ with an upper cut-off E^'^^, and 
summing over the photon polarizations gives the mixing- 
induced CP asymmetry parameter SfCsTr'^'yi^^KTr) 



-2 sin 



in 2/3 1 



nib 



cos ( 



gljUx^) E^ 2 

o7 M\-^ Ktt 

2UgK*K^ E}y ' 



(57) 

n]} 



effect, and the last term is the nonresonant contribution. 
The dependence on the pion energy cut-off is contained 
in 



I{MKrr,Er'') 



I' 



dE^ 



;U2 



(58) 



We show in Fig. 4 results for the Skstt^i parameter as 
a function of Mk-k , integrated with an upper pion energy 
cutoff E^^'^^ = 0.5 GeV. We used in this computation the 
central value for sin 2(3 = 0.685±0.032 as measured in the 
charmonium system The effect of the nonresonant 
contribution is to introduce a mild dependence of the 
asymmetry on Mktt- 

Finally, we integrate also over Mktt = {mK + ti^^, 1.8) 
GeV to obtain the inclusive CP asymmetry parameter 
(for an upper pion energy cut) 



{TTl 
— + hs COS <j)s (59) 

^ Re[KhiEri]} 



The first two terms represent the resonant B K*^ where the phase space factor l2{E™^^) is defined as 



/ dMi^BWj,. (Mk.)|^ f'^^dE^^ 

j^(^max) ^ ^ — = 0.20 + 0.11Z GeV, (60) 

J dMi, JEW K>{Mi,^Wf- dE^\p^\^ 



and the numerical value corresponds to E^^^^ — 0.5 GeV. 
As mentioned, for pion energies above 1.0 — 1.5 GeV, the 
right-handed photon amplitude is power suppressed, so 
it can be expected to be numerically small. We esti- 
mate the contribution from the hard-soft region 0.5 GeV 
< Et; < 1.0 — 1.5 GeV by assuming the validity of the 
low energy expression for the decay amplitudes over this 
range. Taking E™^^ = 1.5 GeV replaces the numerical 
value in Eq. (EOJl with I2 = 0.14-1- O.OSi GeV. In both 
cases discussed above, the contribution of the nonreso- 
nant (third) term in the braces in Eq. (|59|l is less than 
0.5%, and is thus neghgible. 

We neglected in this computation the presence of 
higher K* resonances. The Particle Data Book lists 
four K* resonances in the region niK + 1^1^^ < Mk-k < 1-8 
GeV, which can appear in the Kstt^ invariant mass spec- 
trum (with quantum numbers = 1^, 2+, 3^, . . . ). 
Their inclusion does not change the leading order non- 
resonant right-handed photon amplitude computed here, 



but introduces additional power suppressed effects sim- 
ilar to those parameterized by {hs,4's)- In the narrow 
width limit, their effect is to replace hs cos (ps in Eq. 1(59(1 
with 

hsCOS(j)s (/is cos 0s -fy^Xi /is cos (/)!,) (61) 

with X, = Br{B K*-f)Br{K* ^ KTi)/Br{B K*j), 
and {hi, new parameters for B K*^ defined anal- 
ogously to Eq. 1(43(1 . Of the four kaon resonances con- 
tributing to the sum, only two of them decay into Kt: 
with a branching fraction larger than 30%: ^2(1430) 
and if*(1680), with Xi = 0.15 and 0.01 respectively. 
This shows that the contributions of the higher kaon res- 
onances are likely very small and can be neglected. 

Our results demonstrate that the nonresonant contri- 
bution to the mixing-induced CP asymmetry is negligi- 
bly small, and the SM contamination is dominated by the 
right-handed photon amplitude in B — > K*^fi, parame- 
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terized by hg cos ips ■ Our results show that averaging the 
nonresonant and resonant measurements of the S param- 
eter is a justified procedure. 



research agreement DOE-FC02-94ER40818 and by the 
NSF under grant PHY-9970781 (DP). 



V. CONCLUSION 

We studied in this paper the Standard Model predic- 
tion for the mixing-induced CP asymmetry parameter in 
KsTT^^ decay. This decay is important as a probe 
for new physics manifested through a right-handed pho- 
ton in & ^ 57 decay. The naive expectation for the S 
parameter in the SM is = ~2 sm2(3{ms / nrn) ~ 2%. We 
computed the corrections to this prediction introduced by 
strong interaction effects. 

In the kinematical region with Mk^^^o ^ Mk' and a 
soft pion in the rest frame of the B meson, there is a 
unique SM mechanism contributing to the S parameter 
at leading order in A/ mi,, arising from the B* pole dia- 
grams. These effects are factorizable and calculable using 
a combination of SCET and heavy hadron chiral pertur- 
bation theory In addition, power suppressed effects 
can introduce a potentially sizeable contamination from 
nonfactorizable graphs with the photon coupling to the 
charm quark loop [9j . These are difficult to compute in a 
reliable way, and a simple power counting estimate allows 
a right-handed photon amplitude as large as ~ 9%. 

We performed a detailed numerical study of the non- 
resonant effects. We find that the leading order B* pole 
effect is numerically small. It introduces a weak depen- 
dence of the CP asymmetry SksttOj on Mktt- The domi- 
nant SM contamination is from power suppressed effects 
in the B K*jii resonant amplitude, and our best esti- 
mate in the Mk-r dependent asymmetry is \S^^o^ \ < 8% 
(see Fig. 4). 

When integrated over Mktt, the nonresonant effect is 
practically negligible, and the CP asymmetry SKsTi°-y 
is dominated by the resonant B —f K*jfi amplitude. 
This means that averaging the results of the resonant 
and nonresonant measurements, as currently done at B 
factories, is a justified procedure. If improved measure- 
ments of the CP asymmetry confirm the present average 
\S\ ~ 8%, this would be consistent with a power sup- 
pressed correction in the SM. We reiterate that the naive 
estimate S ~ — 2(r7is/m(,) sin 2/3 seriously underestimates 
the value of the S parameter in the SM. Furthermore, one 
would also expect the agreement between resonant and 
nonresonant measurements to improve. 

We thank Yuval Grossman and Zoltan Ligeti for col- 
laboration on related work and comments. DP thanks 
Jure Zupan for discussions and advice. This work was 
supported in part by the DOE under grant DE-FG03- 
97ER40546 (BG), and by the DOE under cooperative 
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